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ABSTRACT 

Aims. We present a weak gravitational lensing analysis of supergroup SG1120-1202, consisting of four distinct X-ray-luminous 
groups that will merge to form a cluster comparable in mass to Coma at z = 0. These groups lie within a projected separation of 1 to 
4 Mpc and within Av = 550 km and form a unique protocluster to study the matter distribution in a coalescing system. 

Methods. Using high-resolution 7/577ACS imaging, combined with an extensive spectroscopic and imaging data set, we studied the 
weak gravitational distortion of background galaxy images by the matter distribution in the supergroup. We compared the recon¬ 
structed projected density field with the distribution of galaxies and hot X-ray emitting gas in the system and derived halo parameters 
for the individual density peaks. 

Results. We show that the projected mass distribution closely follows the locations of the X-ray peaks and associated brightest group 
galaxies. One of the groups that lies at slightly lower redshift (z ~ 0.35) than the other three groups (z ~ 0.37) is X-ray luminous, 
but is barely detected in the gravitational lensing signal. The other three groups show a significant detection (up to 5(T in mass), 
with velocity dispersions between 35577 q and 530755 km s“* and masses between 0 . 87 q 3 x 10*^' and Lb^gJ x IO'^^/j^'Mq, consistent 
with independent measurements. These groups are associated with peaks in the galaxy and gas density in a relatively straightforward 
manner. Since the groups show no visible signs of interaction, this supports the hypothesis that we observe the groups before they 
merge into a cluster. 

Key words, gravitational lensing: weak - cosmology: dark matter- galaxies: groups: general - galaxies: clusters: general - galaxies: 
formation - X-rays: galaxies: clusters 


1. Introduction 

In the framework of hierarchical structure formation dPeeblesI 
I 1970 I) . matter overdensities grow through merging and accretion 
from the scales of galaxies up to those of large-scale structure 
(LSS). In the concordance ACDM cosmology, the large-scale 
structure of the Universe is driven hy the density fluctuations 
of dark matter, which provide the initial framework for subse¬ 
quent structure formation. As such, the mass distribution in the 
Universe is the driving force behind the formation of cluster¬ 
ing sites for astrophysical processes, such as galaxy groups and 
clusters. 

In galaxy formation and evolution, environment plays a role 
of major imp ortance. Most g alaxies are found in groups and 
clusters (e.g.. lEke et al.l[^04l) . and observations indicate that 
the main part of galaxy evolution takes place in the group en- 
vironment, with signific ant post-processing occurring in clusters 
(iTran et alJl2008i l2009l hereafter T08 and T09). A detailed un¬ 
derstanding of the total mass (dark and visible) and the struc¬ 
ture of the mass density distribution is therefore necessary to 
understand both the processes of group and cluster formation 
and fundamental scaling relations (e.g., iLeauthaud et al.ll2010t 
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iHoekstra et al.l l2012t Ivon der Linden et^ l2014l) as well as to 
distinguish the latter from intrinsic variances in astrophysical 
processes. 


Most overdensities are detected using visible, that is, 
baryonic means. Common methods use g alaxies (red se¬ 
quence and spectrosc opic association, e.g., lEke et al] |2004 
iGladders & Yee 2005I) or gas (X-ray emission or the SZ effect . 


e.g., Sunvaev & Zeldovichl 19701 ri972l:lEinoguenov et al.l2007l) . 

While these methods are efficient, they might not always be as 
effective: they rely on the presence of baryonic matter, while 
the matter distribution is driven by dark matter. Eurthermore, 
the subsequent classification relies on observing the results of 
complex (astrophysical) processes, which introduces a signifi¬ 
cant intrinsic scatter in properties such as X-ray temperatures, 
star formation rate, and galaxy morphologies in different struc¬ 
tures of comparable mass. 

Gravitational lensing is the only direct probe of the total 
mass distribution, in the sense that it does not rely on astrophys¬ 
ical assumptions. A lower signal-to-noise ratio (S/N) makes it 
a less efficient detection method except for massive structures, 
but in combination with complementary methods, it is a pow¬ 
erful independent tool. Erom a statistical perspective, as an in¬ 
dependent, direct measurement, it can s erve as a calibrator for 
mass-observable scaling relations (e.g., ILeauthaud et al.l 120101 : 
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iHoekstra et alJl2Q12l : Ivon der Linden et alJf^l4ll . In individual 
systems it provides an independent estimation of the (projected) 
density field and can shed light on aspects such as interaction, 
dynamical state, and offsets between the overall matter distribu¬ 
tion and the baryonic matter. 

Direct reconstructions of the density fields of individual sys¬ 
tems do not have a high res olution and are predi cted to show 
significant noise fluctuations (Ivan Waerbekdl2000l) . However, it 
gives an important qualitative indication of the dominant density 
distribution, independent of the presence of baryons. To deter¬ 
mine the parameters of the matter distribution, either in a statis¬ 
tical survey or for individual systems, robust centroiding is an 
important task. Using other tracers of the center of mass, such as 
X-ray detections, brightest cluster/group galaxies (BCG/BGG), 
or the luminosity-weighted mean (LWM) position, gravitational 
lensing ca n significantly cons train halo masses and concentra¬ 
tions (e.g.. lGeorge et al.(F2012h . 

Galaxy clusters and, in the past decade, galaxy groups, are 
now identified in a robust manner, including examples of ac¬ 
cretion of smaller structures onto existing clusters. However, 
we have less observational evidence of the connection between 
structures on various scales, that is to say, the initial assembly 
of clusters from groups and galaxies. In t his study, we perform 
a weak lensing analysis of SGI 120-1202 (iGonzalez et al.ll200^ 
hereafter G05), an assembling system of four galaxy groups at 
z ~ 0.37 di scovered in the Las C ampanas Distant Cluster Survey 
(LCDCS, IGonzalez et al.ll200lh . These groups are gravitation¬ 
ally bound and will merge into a galaxy cluster comparable in 
mass to Coma by z = 0. The supergroup, hereafter SGI 120, is 
confirmed by X-ray imaging and optical spectroscopy and has 
already formed a red sequence (see, e.g., G05, T08, T0 9, and 
iKautsch et alj|200^ I Just k al]|201 ll : iFreeland et alj|20flh . 

The individual subgroups are in the low-mass regime of X- 
ray groups, M 200 ~ 10^^ to 10'"^M© and cTy ~ 400 km s“', and 
have not yet interacted. The aim of this study is to determine the 
total matter distribution in the system (dark and baryonic) and to 
constrain individual halo masses. 

This paper is organized as follows. We summarize the gen¬ 
eral framework for weak lensing in Sect. 2. In Sect. 3 we briefly 
describe the data we use, while Sect. 4 covers the framework of 
measurement and analysis methods. In Sect. 5 we discuss the re¬ 
sults and the scientific implications. Section 6 gives a summary 
of our conclusions. 

Throughout this pa per we assume a Planck 
( jPlanck Collaboration et al.l 120141) cosmology with 
Qm = 0.3183, Qa = 0.6817 and Hq - 67.04 km s“^ 
Mpc-'. 


Since we do not know the intrinsic source sizes or magni¬ 
tudes, we can only measure the net distortion or reduced shear 
(^1,^2) = (ri,r2)/(i - /f): 
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where the transformation is written as a multiplication of (1 - k), 
which we do not identify, and a distortion matrix describing the 
alignment of lensed sources in the foreground potential. 

The observed shape of a background source is not a pure 
tracer of the gravitational lensing effect, but the combined ef¬ 
fect of an intrinsic galaxy shape and any distortion of that shape, 
including gravitational lensing. Systematic effects such as tele¬ 
scope aberrations and detector systematics likewise contribute, 
which need to be corrected for. 

The uncertainties the intrinsic shapes introduce is referred 
to as shape noise, and the amount of background sources avail¬ 
able determines the precision of the results. The shape noise 
averages out statistically if we assume that the background 
sources are randomly oriented intrinsically. If faint members of 
the foreground structure are not identified and removed from 
the sample of background sources, these shapes will not be 
affected by gravitational lensing, but might be aligned with 
the potential of the structure under investig ation. These effects 
are k nown as intrinsic alignment (see, e.g.. iMandelbaum et al.l 
l2006h . However, recent results suggest that intrinsic alignments 
should have negligible infl uence for current cluster weak lensing 
studies (ISifon et alJl20L5h . 

The average measured distortion, corrected for systematic 
effects, can then be related to the projected density distribution 
in the lensing structure through 
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where 0 represents the angular coordinates on the plane of the 
sky, 2(0) is the projected density distribution, and are the an¬ 
gular diameter distances between the observer, lens, and back¬ 
ground sources (luminosity distances, sometimes written as D\, 
are not used throughout this paper). 

Normalized by 47rGc“^, the convergence k is therefore also 
known as the dimensionless surface mass density, directly re¬ 
lated to the lensing density distribution and the lensing geome¬ 
try. For axisymmetric lenses, \y\{0) - k{< 0) - k{0) with k{< 0) 
the mean surface mass density within a radial separation 6 - |0| 
to the lens centroid. 


2. Weak lensing framework 

Gravitational lensing is the effect of curved space-time on the 
paths of light rays from distant sources to the observer as they 
pass through the potential of foreground structures. This geo¬ 
metrical effect leads to a displacement of point sources on the 
projected plane of the sky. The differential effect on extended 
sources leads to magnification and distortion effects. This is 
commonly described as a coordinate transformation 
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where the trace component k is known as the convergence and 
the reduced symmetric part is determined by the gravitational 
shear ( 71 , 72 )- 


3. Data 

For this project we made use of results of an extensive multi¬ 
wavelength data set (see, e.g., G05, T08, and T09). 

Key to our lensing analysis are optical data, consisting 
of high-resolution 7/577 ACSQ imaging used for shape mea - 
surements, a s well as VLTyV IMOS (iLe Fevre et alJ l2003l) . 
yL77FORS2 (lAppenzeller et al.l Il998l) . and MagellanfLDSS3 
spectroscopy. 

We also used the X-ray temperatures based on 
ChandralAClS imaging and stellar masses inferred from 

’ Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc. 
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VLTyVIMOS BVR photometry (T08) and complement our opti¬ 
cal color information with ^TPA^O/FLAMINGOS near-infrared 
(NIR) Ks imaging. 

Weusea= 1 i'll 9^5 8^0, d = 12°03'33"0 as center of coordi¬ 
nates, which roughly is the center of the VIMOS imaging data. 
To convert angular to physical separations, we use a reference 
redshift of z = 0.37; this is the median of the redshifts of the 
four BGGs. 



R.A. [deg] 

Fig. 1. Layout of the VLTyVIMOS pointings (red) and 7/577ACS 
pointings (blue). The detected X-ray peaks are shown as well 
(gray), with the radius of the circles 0.5 /t^' Mpc. The X-ray 
peaks 1 and 6 (light gray) are associated with structures at higher 
redshift, beyond SGI 120 (G05). 


3. 1. HST imaging 

The HST/ACS imaging data were taken in July 2005 and January 
2006 and consist of ten pointings, forming a contiguous 11' x 18' 
mosaic. Each tile was observed in F814W (O'.'05/pixel) for 2ks 
over four dithered exposures. Figure [T] shows the layout of the 
different pointings. 

We r educed the data with the same pipeline as em- 
ployed in Schrab back et al.l (l2010l) . which uses MultiDrizzle 
(iKoekemoer et al.l 2003 ) to stack exposures and remove cosmic 
rays. It also includes careful refinement of shifts and rotations 
bet ween exposures as well a s optimized weighting. 

ISchrabback et alJ (l2010h found that using MultiDrizzle 
with the default cosmic-ray rejection parameters can cause cen¬ 
tral stellar pixels to be flagged as cosmic rays, especially when 
there are significant PSF variations between exposures. Galaxies 
are not affected, due to their shallower light profiles. To avoid 
differences in the effective stacked PSFs, we created separate 
stacks for stars and galaxies, with less aggressive cosmic-ray re¬ 
jection for the former. 

For a more detailed descript ion of the reduction process, we 
refer to ISchrabback et al] (|2010|) . 


3.2. Spectroscopy 

We employed optical spectroscopy consisting of three subsets of 
data. The first subset of targets was selected from a magnitude- 
limited catalog (R < 22.5), with preference given to objects in 


visually overdense regions (G05), and observed using VIMOS. 
Follow-up spectroscopy was selected from Ks catalogs (K^ < 
20) and carried out on LDSS3 and FORS2. Figure |3 shows the 
redshift distribution of the final target selection, using a redshift 
quality cut as defined in T08. 



Fig. 2. Redshift distribution of the spectroscopic targets around 
SGI 120. Three significant peaks are identified between 0.35 < 
z < 0.37, associated with SGI 120; between 0.43 < z < 0.44, 
unassociated but concentrated slightly north of peak 1; and be¬ 
tween 0.47 < z < 0.49, associated with X-ray peaks 1 and 6. 


Members of the subgroups of SGI 120 were initially selected 
as galaxies within 500 kpc of their respective X-ray peaks within 
the redshift range 0.32 < z < 0.39. We narrowed the redshift 
range to 0.34 < z < 0.38 without the loss of any members. 
Figure |3] shows a layout of the targets. 




X [h ' Mpc] 

- 3 - 2-10 1 2 3 



-10 -50 5 10 

X [arcmin] 


Fig. 3. Layout of spectroscopic targets, overlaid with the 
VLI/VIMOS pointings. (x,y) - (0,0) corresponds to a = 
ll';i9?'58!0, 6 = 12°03'33'.'0. Colors correspond to the peaks 
in Fig. 13 The BGGs are indicated by larger circles. 
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3.3. Subgroups 


with density profiles to the subgroups based on the BGG and 
X-ray peak positions. 


Table 1. Properties of galaxy groups in SGI 120 and the two 
structures identified at higher redshift. 


BGG 

X-ray 

z 

T 

cr. 

N 

ID 

peak 


(keV) 

(km s“') 


1 

2 

0.3522 

9 9+0.7 
^•^-0.4 

303 ± 60 

13 

2 

3 

0.3707 


406 ± 83 

19 

3 

4 

0.3712 

i-s-i." 

580 ± 100 

29 

4 

5 

0.3694 

-1.0 

567 ± 119 

21 


1 

0.4794 

9 9+0.4 
^•*^-0.3 

820 ± 101 

19 


6 

0.4801 
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In Table [T] we give an overview of the properties of each 
subgroup, also given in G05, T08 and T09. We use the same 
numbering convention for the X-ray peaks and BGGs as these 
papers. Subtle differences in the number of group members are 
due to using the same selection criteria as T09, but a slightly 
different cosmology. 

Figure|4]shows the radial (z, y) distribution of galaxies within 
500 kpc of each subgroup, with 0.315 < z < 0.415. The sub¬ 
group associated with X-ray peak 2 and BGG 1 is at a slightly 
lower redshift, with an estimated 30 ~ 40 Mpc in angular diam¬ 
eter distance to the median redshift of the supergroup. 


Da [h ' Mpc] 

1000 1050 1100 1150 



0.32 0.34 0.36 0.38 0.40 


Fig. 4. Radial (z,y) distribution of all objects (black dots, with 
a redshift quality of 3 as dehned in T08). Objects within 500 
kpc of an X-ray peak and the corresponding selection criteria 
are shown in different colors for easy distinction. The BGGs are 
indicated by larger circles, group redshifts by dashed lines. 


4. Analysis 

In this section we briefly describe our method of shape measure¬ 
ment. We discuss the redshift distribution and selection of back¬ 
ground sources. After establishing a reliable background catalog 
with robust shapes, we describe how we obtain a qualitative re¬ 
construction of the projected mass density and complement this 


4.1. KSB+ shape measurements 


The art of measuring accurate galaxy shapes is an ongoing 
held of investigation, as witnessed, for instance, by the Shear 
TEsting Programmes and the GRavitational lEnsing Accuracy 
Testin g dHeymans et al.l IT oOhl: Massev et al.l 2007I: iBridle et al.l 


testin g dhteymans et al.lUUUot IMassev et al.M/UU/t la ridle et al.l 
I 2 OIOI: iKitching et~ I2012I: iMandelbaum et akl 2015I hereafter 


STEP, STEP2, and GREAT0 8, GREATIO, and GREAT3). We 
make use of the KSB method dKaiser et al.ll995[) . the most com¬ 
monly used and tested technique in the past decade, and discuss 
its application to ACS data. 

Eor this st u dy w e use d the same approach as 
ISchrabback etaP OOO^ IMoi the TS pipelin e in STEP 
and S TEP2) based on the implementation by lErben et ^ 
dMIb- ksb uses the hrst-order effects of distortions induced 
by gravitational shear and PSP on the weighted second moments 
of the light distribution of a source to estimate the reduced 
shear. We detect object s in the same way using SExtractor 
(iBertin & Arnoutslll996h . 

We dehned individual weights for this method based on the 
variance of the shear estimators from this pipeline as 


W cx (o-2„ H- cr2 


+ (T" 



(4) 


where we assumed a minimum variance of cr^„, the intrinsic 
shape noise. Based on the hndings in the STEP analyses, we 
expect an underestimation of the shears by KSB-t of about a 
few percent. We applied the same empirical correction factor as 
ISchrabback et al.l (120lOh to account for this expected bias. 

The systematic distortion effects due to telescope and optic 
system give rise to shapes convoluted by a point spread function 
(PSF). The main source of variations of the ACS PSF is caused 
by changes in the telesco pe focus, causing spatial and tempo¬ 
ral hu ctuations (see, e.g.. ISchrabback et alJl200'^ iRhodes et al.l 
I2nn7h . 

A common strategy is to map the distortions caused by the 
PSF using the shapes of foreground stars, but the average num¬ 
ber of stars in our ACS images is ~ 20 - 40. This leads to a 
poorly sampled PSF and an imperfect correction, causing sig- 
nihcant residual distortions, especially detectable at the edges of 
the images where the tiles overlap sligh t ly. We therefore adopted 
the same strategy as ISchrabback et al.l (1201 Oh based on a prin¬ 
cipal component analysis of the PSF variation in dense stellar 
helds. 

Furthermore, deterioration of the ACS CCDs over time due 
to constant exposure to cosmic rays in space leads to an ef¬ 
fect called charge-transfer inefficiency (C TI), causing charge 
trails in the CCD readout direction (e.g., [Rhodes et all 120071 : 
iMassev et al.ll2010t ISchrabback et alJl2010h . These effects will 
affect the measured shear pattern and the reconstruction of the 
projected density distribution, and it is therefore important to 
correct for them. 

Here we applied the same par ametric CTI model as de¬ 
scribed in ISchrabback et alJ (l2010h for the correction of the 
KSB-i- polarizations. 

After constructing the shape measurement catalogs, we ap¬ 
plied several common selection criteria and cuts. These crite¬ 
ria are based on simulations and quality flags of the detection 
and shape measurement pipelines, and they depend on the noise 
properties, on the variance and convergence of the model hts, 
and on the object and PSF size. 
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A list of the various selection criteria can be found in the 
appendix. Sources that pass the criteria of all pipelines num¬ 
ber 7012, for a source density of ~ 64 galaxies/arcmin^ with 
MAG-AUTO magnitudes tgu < 27.1. 


4.2. Redshift distribution 

We acquired spectroscopic redshifts for 497 objects in our cat¬ 
alogs. The spectroscopic targets were selected based on mag¬ 
nitude, and preference was given to visually overdense regions, 
which means that these spectroscopically conhrmed members do 
not give a complete picture of the galaxy distribution in SGI 120. 
The brightest conhrmed supergroup member has /gw = 17.5, 
while the spectroscopic survey remains > 50% complete to 
/gw = 20.5 (T09). We hnd that conhrmed supergroup members 
have numbers peaking between magnitudes 19.5 < /gw < 20.0. 

To separate background and foreground sources, we con¬ 
sidered that group members are expected to dominate number 
counts in the magnitude range of conhrmed members, while 
background number counts dominate at fainter magnitudes. We 
initially selected background sources as objects with /gw > 22 
and assessed possible contamination by faint foreground objects. 
Figure|5]shows the number density of sources with /gw < 22 and 
/gw > 22, where we have used a Gaussian smoothing with a 20” 
kernel width. 

Because gravitational lensing is a geometric effect that has a 
non-linear dependence on redshift, we took the expected redshift 
di stributions into account, f ollowing the same parametrization as 
in ISchrabback et all (l2010l) . We show the total redshift distribu¬ 
tion of sources with /gw > 22 in Fig.|2 For a given lens redshift, 
such as in this particular system, the lensing signal has a lin¬ 
ear dependence on the lensing efhciency = max {0, Z)is/Dos}. 
We can therefore determine a mean lensing efficiency {y6) for the 
sources with respect to each subgroup redshift. 

As mentioned earlier, both X-ray peak 1 and 6 (G05) are as¬ 
sociated with structures at higher redshift (both 0.46 ^ z ^ 0.48). 
We must take the gravitational distortions caused by these back¬ 
ground structures into account when trying to isolate the signal 
from SGI 120. We therefore also determined a mean lensing ef¬ 
hciency for these structures. 

We found average lensing efficiencies of (^6) » 0.52 for 
SGI 120, corresponding to an effective background redshift of 
Zeff ~ 0.88, and (/?) s; 0.42 for the two background structures, 
corresponding to an effective background redshift of Zeft ~ 0.95. 


4.2.1. Foreground contamination 

An intrinsic redshift distribution of sources with /gw > 22 im¬ 
plies that some of these objects are faint foreground sources 
or members of the SGI 120 structure. Based on our parametric 
redshift distribution, we estimate that ~ 9% of our background 
sources to lie in front of SGI 120. 

Foreground sources are not lensed by the groups. We account 
for this dilution effect by assigning y6 = 0 to this part of the 
redshift distribution in our dehnition of the lensing efficiency 
above. 

This assumes a random held of view, which is not the case 
for our observations, with known overdensities at z ~ 0.37 and 
z ~ 0.48. However, Fig. |5] suggests no signihcant correlation 
between the distribution of these sources and the galaxy distri¬ 
bution of SGI 120. To estimate possible variations in the number 
density n of sources with /gw > 22, we derived an average num- 



0.0 0.2 0.4 0.6 0.8 1.0 


Fig. 6. Parametric redshift distribution of sources with /gw > 22 
(upper panel) and the corresponding distribution in lensing ef¬ 
hciency p (lower panel). In the upper panel, the dashed green 
line corresponds to Zeft ~ 0.88 with respect to SGI 120 and the 
dotted red line corresponds to Zeft ~ 0.95 with respect to the 
two structures at higher redshift. In the lower panel, the dashed 
green curve shows the distribution in p with respect to SGI 120, 
with (P) Si 0.52 (dashed green vertical line), and the dotted red 
curve the distribution in p for the two structures at z = 0.48, with 
(P) s! 0.42 (dotted red vertical line). 


her density prohle around the group centers, as shown in Fig. 

m 



Radial separation [arasec] 

Fig. 7. Variations in galaxy number density n as a function of ra¬ 
dial distance from the lens positions, using the X-ray peaks (di¬ 
amonds) and BGG positions (squares). Data points are slightly 
offset for clarity. Overplotted is the average number density of 
~ 64 galaxies/arcmin^ of the whole ACS mosaic (black solid 
line) and the best-ht radial prohle (dashed) with Icr errors (dot¬ 
ted). The estimated effect of the lensing magnihcation, is 
shown in grayscale, varying the slope of the luminosity function 
between 0 < cr < 3. Different shades in grayscale correspond to 
steps of 0.5 X 10'^ in group mass M2oo- 


We used radial bins between 10” < 0 < 95” to avoid the 
BGGs and the edges of the ACS coverage. We also considered 
only the group centers of groups 2 through 5, as group 6 is not 
entirely covered by the ACS pointings. Finally, we averaged the 
signal over all four subgroups to increase the signal-to-noise ra¬ 
tio. 

We then quantihed the radial dependence of the galaxy num¬ 
ber density by htting a parameterized prohle given by n(6) - 
(1 -t al6)nbg, with 0 in arcseconds and with Wbg = 64/arcmin^ 
hxed. (In fact, if we allow Ubg to vary, we recover Wbg = (64 + 
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3)/arcmin^.) We found a = 0.14 + 1.11, consistent with no trend 
in galaxy number density with radial separation from the lens 
centers. 

To interpret this radial number density profile, we have to 
consider both the presence of unidentified faint group mem- 
bers and the effect of the lens ing magnification fi (see, e.g., 
iBartelmann & Schneider! 1200 iL and references therein). The 
presence of unidentified group members would increase the 
number density. Magnification increases both the observed flux 
of background sources, leading to an increase in n, and the solid 
angle behind the lenses, causing a dilution of n (not to be con¬ 
fused with the dilution of the shear signal caused by unlensed 
foreground objects in the background source sample). It then 
depends on the slope a of the luminosity function whether the 
lensing magnification causes a net increase or decrease in num¬ 
ber density by *, where fi and a depend of the source r ed- 
shift. Both effects were shown by iHildebrandt et al.l (l2009l) . A 
decrease could cancel the effect of unidentified group members. 

We wish to obtain a rough estimate of the expected influence 
of magnification to check whether it is smaller than the statis¬ 
tical uncertainty. For this we ignored the redshift dependencies 
of yu and a and considered a wide range 0 < a < 3, which was 
simply chosen to assess all possible variations in the magnifica¬ 
tion without making assumptions about the luminosity function. 
We used a group mass of M 200 = 1.0 x lO'"^, where M 200 is de¬ 
fined as the total mass within a radius of r 2 oo of a halo, where the 
mean density of the halo is p(< r 2 oo) = 200pcrit(z). For a robust 
estimate, we also considered variations in M 200 of +0.5 x lO''^. 
The results are plotted in Fig.|2]in grayscale. 


4.3.2. Density profile parameters 

Earlier studies indicate that the groups are infalling for the first 
time and have not yet interacted, although X-ray measurements 
show a possible onset of interaction (G05). We considered the 
groups as individual overdensities with spherically symmetric 
density distributions and derived halo parameters for each group, 
including the background structure around X-ray peak 6. 

We considered two types of density profiles and two possible 
choices of gr oup centroids. We c onsidered the Navarro-Frenk- 
White (NFW, iNavarro et al.lll99^ density profile and compared 
this to the singular isothermal sphere (SIS) model. 

The SIS profile is determined by a single free parameter, the 
halo velocity dispersion cTy, where the subscript y is used to dis¬ 
tinguish this parameter, derived from a two-dimensional model 
of the projected mass density, from other derivations of velocity 
dispersion, such as the one-dimensional ctj. derived from the red- 
shift distribution. The advantages of this profile are its simplicity 
and the linear dependence of the lensing signal on the squared 
velocity dispersion. The tangential component of the shear with 
respect to the group center is given by 


rt(fl) = 


In 2/3 


(5) 


where 6 indicates the separation from the center in radians. This 
allows for a straightforward interpretation of any possible corre¬ 
lation between the fitted parameters of different subgroups. 

The NFW profile is usually expressed in terms of its mass 
and concentration and depends on redshift. The halo mass M 200 
is given by 


We estimate possible magnification effects to be smaller than 
the statistical uncertainties from potential residual group mem¬ 
ber contamination. We expect that any residual excess contami¬ 
nation by member galaxies of SGI 120 in the source sample must 
be small and comparable to regular line-of-sight variations. In 
Sect. 5.2, we confirm that we do not need to apply a dilution 
correction for excess contamination from the supergroup itself. 

4.3. Lensing analysis 

Our approach to determining the matter distribution in SGI 120 
is twofold. First, we show that the distribution of light (galaxy 
number densities, BGGs, and X-ray peaks) is closely correlated 
with the underlying mass distribution. Second, we determine the 
density profile parameters for each subgroup, taking into account 
the effect of each subgroup and background structure simultane¬ 
ously. 

4.3.1. Reconstruction of the mass distribution 

We used a Kaiser-Squires (KS, iKaiser & SauiresI flOOSh in¬ 
version technique to reconstruct the surface mass density. We 
smoothed the data onto a rectangular grid, using a Gaussian 
smoothing kernel with a width of 20”, equal to the smoothing 
used for the galaxy number densities in Fig. |5] 

We investigated possible systematic errors in our data by 
changing the phase of the shear by jtt, which corresponds to 
rotating the background galaxies by ^tt. The distortion caused 
by weak lensing does not introduce a curl in the shear held, and 
the resulting reconstructed map should display only noise in the 
absence of systematic errors. 


M 200 = 200pcrit(z)^^''2oo = 100^^r2oo- (6) 

The concentration C 200 is defined as the relation between the 
characteristic shape of the density profile and r 2 oo- The analyti- 
ca l formulas for the shear sig nal of an NFW profile c an be found 
in lWright & BrainerdI(l2000li and iBartelmannI (Il996l) . 

Because of the lower S/N, the centers of dark matter haloes 
should not be estimated directly from the lensing data when de¬ 
termining density profile parameters. Instead, one has to rely 
upon visible tracers such as peaks in the X-ray emission of hot 
gas or the brightest or heaviest galaxy (e.g., in terms of a stel¬ 
lar mass as derived in T08) in the group or cluster. If the fitted 
halo model is offset from the true underlying halo, the fit is in¬ 
ferior and the introduced systematic uncertainties can be signif¬ 
icant (iGeorge et al.ll201^ . In particular, the halo mass will on 
average be underestimated, while the uncertainties, most often 
determined from confidence levels, will be increased. This leads 
to both a biased and a less effective st udy. 

As described in iGeorge et al.l (l2012h . there are several 
choices possible as tracer of the halo center. These can be based 
upon a central galaxy, several or all of the associated galaxies, or 
the X-ray flux. In this study, the haloes under consideration are 
part of a coalescing system, and an offset from the true halo cen¬ 
ter of some or all of these tracers is not unlikely. However, the 
BGGs of the subgroup are als o the most massive group galaxies 
(MMGG, IGeorge et al.ll2012l) in terms of stellar mass and mag¬ 
nitudes in most observed bands and coincide well with the X-ray 
peaks (T08). We derived the parameter values using both options 
and determined whether these are consistent. 

Given the close angular separation of the X-ray peaks, we 
did not compute azimuthally averaged profiles. Instead, we com¬ 
puted the total lensing distortion g - Yi Si for sach background 
source induced by each of the six foreground structures. This is 
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valid if we assume g « 1, which is certainly the case for the 
sources where the distortion is not dominated by one of the lens- 
ing structures. 

We then determined prohle parameters for each subgroup us¬ 
ing minimization. For X-ray peak 1, we assumed cr - 820 
km s“^ from G05 and an order of magnitude M 200 = 3.7 x 
and assessed the effect of omitting the influence of 
this background structure. 


5. Results 

In this section we discuss the reconstructed density distribution 
and best-ht prohle parameters, and we show that SGI 120 is con¬ 
sistent with expectations from hierarchical structure formation, 
even though the system is not relaxed. 

5. 1. Matter distribution 

In Fig. 0 we show the reconstruction of the projected surface 
mass density. We detect signihcant peaks near three of the fore¬ 
ground structures. We do not detect a signihcant peak in the den¬ 
sity distribution near X-ray peak 2. 

We considered the results of our mass reconstruction in a 
qualitative manner. The peaks in our surface mass density re¬ 
construction coincide very well with the peaks in galaxy num¬ 
ber density (Fig. |5]l and X-ray emission (G05), within smooth¬ 
ing scales. We found no signihcant ‘dark’ overdensities either, 
and small offsets between peaks using various tracers are ex¬ 
pected in a coalescing system. Finally, the map shows signih- 
cantly stronger peaks than the control map. 

5.2. Individuai groups 

5.2.1. SIS velocity dispersions 

We present the results of the joint minimization ht of SIS 
prohle parameters around the X-ray peaks in Fig.|9] The reduced 
X^ value isxl - ^ .4. 

The combined contours of Fig. |9] show no features that in¬ 
dicate signihcant degeneracies between the individual group cr^ 
values. While it is to be expected that nearby mass concentra¬ 
tions inhuence the shear pattern around an individual lens, we 
conclude that noise is a dominant factor in these results. More 
massive haloes or smaller halo separations can be expected to 
increase correlations. 

The resulting cry values are given in Table|2l Consistent with 
the reconstructed mass map in Fig. 0 we do not detect a very 
signihcant lensing signal around X-ray peak 2, barely exceeding 
the 68% conhdence limit. 

The velocity dispersion associated with X-ray peak 1 is nec¬ 
essarily kept constant, as the peak lies outside the ACS mosaic. 
Upon inspection, it turns out that varying this parameter between 
0 < (Ti < 820 km s“* does not alter the results by more than 10% 
of the 68% conhdence interval for X-ray peak 2, which lies clos¬ 
est to peak 1. The effect is even smaller for the other groups. 

Similar to our assessment of systematics for the mass map 
reconstruction, we repeated the ht to a control signal by chang¬ 
ing the phase of the shear by jtt. The results are consistent with 
a control signal of » 0. Because of its less favorable lensing 
geometry ((J3) - 0.42), the constraints for group 6 are weaker, 
although it is still detected at a signihcance of cr » 1.6. 

Finally, we determined how much our results would be af¬ 
fected if the signal were boosted by a dilution factor of 1 -H (a -H 


Peak 2 



[100^ km^ s“^] 


Fig. 9. Marginalized 2D x^ distributions of the simultaneous ht 
to the individual subgroup velocity dispersions, together with the 
marginalized ID likelihoods for each subgroup. Overplotted are 
the 68.3%,95.4%, and 99.7% conhdence levels. 


CTflj/fl = 1-1-1.25/6 for group member contamination, as dis¬ 
cussed in Sect. 4.2, using a conservative Icr upper limit. We hnd 
that this does not alter the results by more than 37% of the 68% 
conhdence intervals, justifying our earlier approach. 

We repeated the ht around the BGGs as tracers of the halo 
centers. The results are very similar, with the htted values also 
given in Table |2] There is some difference with up to 2cr devia¬ 
tions between the results for peaks 3 and 5, where the separation 
between X-ray peak and BGG is also the largest. The quality of 
the ht, in terms of a reducedvalue, is the same. 

5.2.2. M 200 

In the same manner, we determined NFW prohle parameters 
from the distortion pattern in the ACS held around the sub¬ 
groups. 

Weak lensing data of individual groups or low-mass clus¬ 
ters do not have sufficient signal-to-noise to provide useful 
constraints on M 200 and C 200 simultaneously. Therefore, 
we employed th e mas s-concentration relation given in 
iMandelbaum et all (l2008l) . restricting the ht to one free 
parameter, M 2 oo- The results of these hts are summarized in 
Table |2] both for the X-ray centroids and BGGs as tracers of the 
halo centers. 

5.2.3. Scaling relations 

G05 showed that the subgroups w ere consistent with the local 
Tx - cr^ relation (IXue & Wull2000l) . a fact which did not change 
with more spectroscopic data in T09. Here we did not deter¬ 
mine ID velocity dispersions from the redshift distribution of 
group members, but assumed the projections of 3D halo models. 
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Table 2. Profile parameter fit results 


Subgroup 

O'z 

cTy (X-ray) 

ay (BGG) 

4^200 (X-ray) 

4^200 (BGG) 

(X-ray ID) 

(km s“*) 

(km s“*) 

(km s“') 

(lO'Vi-'Mo) 

(lO'Vi-'Mo) 

2 

303 ± 60 

240+90 

230+95 

u.D_o.2 

^•^- 0.2 

3 

406 ± 83 

530:^5 

49^4-60 

^•^-0.4 

^•^- 0.3 

4 

580 ±100 

-50 

445!” 


-^•^- 0.4 

5 

567 ± 119 

355)” 

480!” 

0 . 8 !“:^ 

^•^- 0.4 

6 


325^^^ 

305!^°5 

o.7!°:S 

0 


Hence, we are not limited by group member identification. As 
mentioned before, group centroiding can be a problem. 

Although the parameters of individual groups have shifted in 
this analysis, on average the groups still lie on the local Tx - cr^ 
relation, showing a scatter of similar magnitude as the data in 
IXue & Wul (120001 Fig. [ig. 




Fig. 10. Comp arison of the pr operties of SGI 120 with ob- 
served <t - T (|3^e & WulbOOOl left) and Lx - M 200 relations 
dLeauthaud et ^ 2010 . right). Results based upon X-ray peaks 
and BGGs as center of mass are indicated by crosses and open 
squares, respectively, while velocity dispersions from G05 and 
T09 are shown as open diamonds. Horizontal error bars are plot¬ 
ted at the vertical median. 


iLeauthaud et al.l (1201 Oh constrained the Lx - M 200 scaling 
relation using weak lensing data of groups in the COSMOS field. 
The supergroup as a whole is consistent with this scaling relation 
as well, within the scatter (Fig. [ig. 

These interpretations would be reinforced if we did not take 
X-ray peak 2 into account and considered the conclusion in 
iGeorge et al.l (l2012h that BGGs/MMGGs are better tracers of 
group halo centers than X-ray centroids. 

Even though individual groups do not always lie precisely 
on the determined scaling relations, differences in environment 
and their effect on the astrophysical processes behind the observ¬ 
ables used in these analyses create intrinsic scatter around these 
relations, w hich is averaged out in a stacking analysis such as 
employed in lLeauthaud et alJ(l2010h . 


6. Summary 

We have performed a weak lensing analysis of the coalescing su¬ 
pergroup SGI 120 and showed that the underlying density distri¬ 
bution of matter is well traced by both visual galaxy light and X- 
ray emission. The subgroups of SGI 120 have not yet interacted, 
but are expected to do so within short timescales, as projected 


separations are of about 1-4 Mpc (G05). As such, the system is 
a unique demonstration of hierarchical structure formation. 

Slight offsets between the peaks in the galaxy distributions. 
X-ray gas, and the total matter distributions are well within 
smoothing scales used and are consistent with an unrelaxed sys¬ 
tem on the verge of merging. We found that us ing either X- 
ray p eaks or BGGs as tracers for the halo centers (iGeorge et al.l 
1200 has a minor impact on the derived halo parameters, with 
results consistent within 2cr error bars. We consider these con¬ 
clusions to be an indication of the robustness of our results. 

Furthermore, while the groups are close enough to be gravi¬ 
tationally bound (G05), the individual group halo masses are low 
enough compared to their separations to treat them as individual 
lenses, within parameter error bars. 

The fitted profile parameters are consistent with well- 
demonstrated scaling relations, within the intrinsic sca tter cre¬ 
ated by astrophysical variations (ILeauthaud et alJl201flh . This is 
further conhrmation that the observed structure of SGI 120 is 
consistent with the paradigm of hierarchical structure formation, 
providing a unique example of this theoretical framework. 

Structures such as SGI 120 are rare. In fact, SGI 120 should 
be seen as a single piece of a much larger puzzle, where con¬ 
firmation from studies of similar structures is a necessity. The 
structure of SGI 120 is uniquely oriented in the plane of the 
sky, and the subgroups show no signs of interaction yet, mak¬ 
ing it well suited to distinguish the various components and 
overdensities. An example of a well-studie d heavier structure is 
the Cl 1604 supercluster (iGal et al.ll2008ll . where the complex 
structure presents difficulties in determining accur ate masses, ei¬ 
ther u sing spectroscopic velocity dispersions (e.g.- lLemaux et all 
I 2 OI 2 I) or weak lensing analyses of a fe w selected subclusters 
(iMargoniner et alJl2005l:lLagattutall20Tll) . 

Especially the extension of studies like these to individual 
systems of lower mass like SGI 120 will present a significant 
challenge, both in detecting such rare coalescing systems and 
in obtaining robust and accurate lensing measurements, given 
the lower S/N. An interesting approach is the combination of 
large existing spectroscopic group catalogs (e.g., Eke_eLakl2QQ 
iBerlind et all 120061 : iTenipel et alj|2012l : [Robotham et aklboi ll 
and recent or currently ongoing large sky imagin g surveys of var- 
ious w idth and depth, designed for lensing (e.g.,iHevmans et al.l 
IMl iGilbank et al] IMU Ide Jong et alJ |201^ that are sup¬ 
ported by extensive spectroscopic and color information. 
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Appendix A: Quality flags and selection criteria for 
background sources. 

We assigned several quality flags to the source catalogs during 
detection and shape measurement. 

We use d the same rms noise m odel and deblending pa¬ 
rameters as ISchrabback et alJ (1201 Ol) for object detection with 
SExtrator. In addition to detection flags, we required at least 
eight adjacent pixels with values more than 1.4cr above the back¬ 
ground. We defined an initial S/N cut by flagging objects with 
FLUXJVUT0/FLUXERR_AUT0< 10. 

We furthermore selected sources with a minimum size com¬ 
pared to the smearing induced by the PSF. We excl uded sources 
for w hich the half-light radius r/, (as defined in lErben et ^ 
1200 Ih compared to that of the average star is not smaller than 
r/, > 1.2r*. 

Finally, we selected sources w ith a KSB shape measurement 
S/N (defined in lErben et al.ll200Tl) larger than 4, to be consistent 
with KSB-i- studies using a similar definition of the source S/N. 
In this pipeline, the effect of smearing and shearing by the PSF is 
for an important part described by the tensor. To avoid being 
domina ted by noise, we e xcluded sources for which Tr(P^)/2 < 
0.1 (see lErben et al.l200lL for technical details and terminology). 

In the final source selection, the catalog of 8273 galaxies is 
reduced to 7012, ~ 64 galaxies/arcmin^, that pass all quality cri¬ 
teria from detection and shape measurement. 
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Fig. 8. Smoothed map of the reconstructed projected density distribution (left) and the imaginary control signal (right) where the 
shear signal is rotated out of phase. X-ray peak positions are indicated by white crosses and BGG positions by white circles. 
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